Gravitational wave burst search in the Virgo C7 data 
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A search for gravitational wave burst events has been performed with the Virgo C7 commissioning 
run data that have been acquired in September 2005 over five days. It focused on un-modeled short 
duration signals in the frequency range 150 Hz to 2 kHz. A search aimed at detecting the GW 
emission from the merger and ringdown phases of binary black hole coalescences was also carried 
out. An extensive understanding of the data was required to be able to handle a burst search using 
the output of only one detector. A 90% confidence level upper limit on the number of expected 
events given the Virgo C7 sensitivity curve has been derived as a function of the signal strength, 
for un-modeled gravitational wave search. The sensitivity of the analysis presented is, in terms of 
the root sum square strain amplitude, hrss — 10"^"/^ Hz. This can be interpreted in terms of a 
frequentist upper limit on the rate IZgo% of detectable gravitational wave bursts at the level of 1.1 
events per day at 90% confidence level. From the binary black hole search, we obtained the distance 
reach at 50% and 90% efficiency as a function of the total mass of the final black hole. The maximal 
detection distance for non-spinning high and equal mass black hole binary system obtained by this 
analysis in C7 data is ~ 2.9 ±0.1 Mpc for a detection efficiency of 50% for a binary of total mass 
8OM0. 
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PACS numbers: 04.80.Nn, 07.05Kf 



1. INTRODUCTION 

Virgo [l| is a 3-km long arm power-recycled Michelson interferometer located near Pisa, Italy, whose goal is to detect 
gravitational waves (GW) emitted by astrophysical sources extending out past the Virgo cluster. The commissioning 
of the detector started in 2003 and regular data taking campaigns have been organized after each important milestone. 
The last commissioning run (C7) took place in September 2005 and lasted for 5 days. The best achieved sensitivity 
was h ~ 7 X 10~^^/\Affz at 300 Hz. The Virgo design sensitivity at this frequency is expected to be better by an order 
of magnitude assuming that 10 W enters into the interferometer. However, during the C7 run, Virgo was running 
with a reduced light power, 0.7 W, because the backscattering in the mode-cleaner cavity of the light reflected by 
the recycling mirror prevented the control of the interferometer at full power. Despite the reduced sensitivity, several 
GW searches have been carried out using this data set in order to set up and tune the pipelines developed in Virgo: 
the search of the coalescence of neutron stars the search of continuous wave sources Q, a joint Virgo, AURIGA, 
EXPLORER and NAUTILUS GW search [|, a targeted GW burst search associated with the long GRB 050915a ^, 
and the search of GW burst signals that is reported in this paper. 

The category of gravitational wave bursts includes all possible signals whose duration is short, less than a few 
hundreds of milliseconds. Many violent astrophysical phenomena will be accompanied by an emission of GW burst. 
If these events happen sufficiently close by, the ground based gravitational detectors will be able to observe these 
short duration GW bursts. There are many candidate GW burst sources. These include massive star core collapse 
merging phase of coalescing compact binary systems forming a single black hole (BH) [Til. [T^. [isl [l^ . 
BH ring-down lisll . a strophvsical engines that generate gamma-ray bursts (GRB) [l^, neutron star oscillation modes 
and instabilities [ij], or cosmic string cusps and kinks [Iq. Some of these sources are well modeled, but not all, and 
as such a burst pipeline is built by making very few assumptions about the nature of these waveforms. The theoretical 
event rates for many of these sources are quite uncertain. 

Binary black hole (BBH) mergers are an interesting source for GW burst searches. The typical event rate for these 
sources is 1 Myr~^ per galaxy [H, [1^ and is still highly uncertain given the lack of direct observational evidence for 
the existence of BBH systems, unlike the double neutron star systems. Though these sources are routinely searched for 
by the inspiral phase prior to merger using matched filtering techniques, an independent search focusing on the merger 
and ring-down GW emission, which is dominant for high mass binaries [ll| , may bring additional confidence [63| • A 
search of this type was first attempted with the LIGO S2 data ; this used the then-available numerical relativity 
waveforms from the Lazarus project as a model of BBH mergers. However the estimates obtained were projected as 
order of magnitude estimates due to the nonrobusteness of the waveforms used. Though the burst searches do not 
have as large of a distance reach as an inspiral search, it will be interesting to understand the implications of a BH 
binary search from a burst perspective, especially because a burst search is sensitive to the merger and ring-down 
phase of a BBH coalescence. Among the three phases of inspiral, merger and ring-down, most of the energy may be 
released in the highly relativistic phase of merger, that is difficult to model by analytic approximation methods. But 
recently numerical relativity simulations have made significant advances in generating the waveforms for all the three 
phases of the binary evolution (see [22| and references there in). 

In this paper, we report on an all-sky burst search for un- modeled waveforms and BBH mergers without using any 
prior information on the expected waveform. This concerns all short duration (^ 1 s) signals having energy in the 
best sensitivity frequency band of the C7 data, 150-2000 Hz. This search has been performed using the data of only 
one interferometer which does not allow one to perform timing coincidence to eliminate _^urious events in one of the 
detectors, as was done in the burst searches performed using the LIGO detectors [111, [H, and the bar detectors 



25 1 . However, a similar one-detector burst search has already been carried out with the data from the TAMA detector 
26|. The C7 run data has been used for extensive analysis of the Virgo noise; this is a fundamental step in the path 
toward performing a GW burst search with only one detector. Indeed, numerous sources of noise generate, in the 
interferometer's GW strain amplitude channel, transient events which mimic GW burst events; environmental noise, 
such as acoustic noise, has been found at the origin of many of them. Non stationary data can also generate an excess 
of short duration events. Such phenomena have been found, examined, and understood in the C7 data. These studies 
allowed us to define the pre-processing and post-processing steps needed in order to optimize the performance of the 
GW burst pipeline used in this analysis. 

The outline of this paper is as follows. In Section[2]we describe the main features of the Virgo detector, emphasizing 
the operating characteristics during the C7 run that played a role in understanding the quality of the data. In Section 
12.21 we describe how we select the data periods that have been used for this analysis. The burst pipeline used in this 
analysis is described in Section [3) In Section 13.41 the main results on the C7 data characterization needed in order 
to understand and to suppress the sources of glitches are summarized; this includes the veto strategy used against 



3 



identified sources of noise. Section [4] gives the results of the search in the C7 data. Section [5] explains how the analysis 
sensitivity has been estimated considering different types of possible GW burst waveforms. We especially considered 
in this paper BBH merger and ring-down numerical waveforms [l^]. Section [6] gives the upper limits obtained at 90% 
of confidence level on the number of events as a function of the signal strain amplitude. We finally conclude this 
paper with an astrophysical interpretation of the present GW burst search. 

2. VIRGO DURING THE COMMISSIONING RUN C7 
2.1. Detector status 

The Virgo detector \2E\ is a power recycled Michelson interferometer with 3-km long arms that each contain a 
Fabry-Perot cavity. All mirrors are suspended from the so-called Superattenuator [29'|, whose goal is to reduce 
drastically above 10 Hz the seismic noise transferred to the instrument. A 20 W Nd:YAG laser is used to illuminate 
the interferometer. The laser light is modulated in phase at the frequencies ~ 22 and 6 MHz; this technique permits 
the GW strain to be detected at the modulation frequency where the laser power fluctuation is much smaller than 
in the interferometer bandwidth. The beam is spatially filtered with a 144 m long triangular input mode-cleaner 
cavity before being injected into the main interferometer. The laser frequency is pre-stabilized in order to acquire the 
control of the different optical cavities, but to reach the extreme sensitivity targeted by Virgo an enhanced control 
of the laser frequency noise is required; it has to be reduced by several orders of magnitude. This is the role of the 
so-called second stage frequency stabilization which is engaged during the cavities' lock acquisition ^3d\. The beam 
entering the interferometer is divided by the beam splitter (BS) into two beams that are injected into the 3-km long 
arm cavities. Apart from the mirrors' losses, all light fed-in by the injection system subsequently returns to it. The 
power-recycling (PR) mirror, with a reflectivity of 92%, reflects the out-going light back to the main interferometer. 
Together with the Michelson interferometer the power-recycling mirror forms a Fabry-Perot-like cavity in which the 
light power is resonantly enhanced, thereby improving the shot noise limit. The Michelson interferometer is held on 
the dark fringe, and the GW strain signal is expected in the beam at the dark port, which leaves the vacuum via the 
so-called detection bench. The detection bench is a suspended table accommodating several optical components. The 
beam coming from the BS passes through an output mode-cleaner, a 3.6 cm long rigid cavity. The main output beam 
is detected by a pair of InGaAs photodiodes. Useful signals are obtained by detecting the light in the transmission of 
the arm cavities and in reflection of the power-recycling cavity (B2 photodiode). The GW signal that results from a 
detuning of the carrier resonance in the arms is extracted from the dark port channel, demodulated, and then sampled 
at 20 kHz. This signal is digitized and filtered. 

The control of the interferometer consists in maintaining the laser light resonant in the optical cavities and the 
output port tuned on the dark fringe, defining its working point. More precisely, the carrier must be resonant in all 
cavities while the sidebands must be resonant in the central cavity but anti-resonant in the arms. Despite the good 
seismic noise attenuation provided by the Superattenuator, feedback controls are mandatory in order to keep the 
interferometer locked on the right working point. In addition to the control of the longitudinal degrees of freedom of 
the cavities, the mirrors must be kept aligned with respect to each other. This is the so-called auto-alignment control 
(AA). If uncontrolled, the angular degrees of freedom of the suspended optics distort the cavity eigen-modes. This 
causes power modulation of the light fields; furthermore, long term drifts will make the longitudinal control impossible 
after a certain amount of time, and misalignments increase the coupling of other noise sources into the dark port. 
This has been a major problem in the C7 data analysis (see Section [H^ . In actuality, during C7 not all of the mirrors 
were under AA control; the arm input mirrors and the injection bench were only controlled locally. 

To convert the signal received on the dark port into a measure of GW strain one needs to calibrate the residual 
arm-length difference; this includes at low frequencies the effect of the control loop keeping the cavities in resonance. 
The loop correction signals are subtracted from the dark port signal that is then converted into a strain via knowledge 
of the optical gain of the feedback control loops. This gain may vary depending, for instance, on the alignment drifts 
of the optical elements. To monitor the variation, some sinusoidal signals (four between 100 and 110 Hz, and four 
between 350 and 360 Hz) are applied on the end-arm mirrors and the amplitude of the lines are measured. Another 
important aspect of the GW strain reconstruction is the good knowledge of the transfer function of the actuator chain 
which is used to control the mirrors. Actually, during the C7 run, the calibration procedure was not as accurate as 
it is now. For instance, the frequency dependence of the digital to analog converter transfer function was not taken 
into account in the reconstruction of the GW strain amplitude channel h(t). This yielded a systematic error on h{t) 
which has been estimated to be 40% [31] . 

Moreover, during C7 Virgo was running with a 700 mW laser beam entering the interferometer in order to avoid 
instabilities due to the backscattering of light into the mode-cleaner cavity reflected by the recycling mirror. The 
problem has since been fixed during a shutdown carried out just after the C7 run. This reduced power light limited 
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FIG. 1: Sensitivity curves obtained between 2003 and 2005 during the commissioning of the Virgo detector before the 2005 
shutdown. The black curve is the Virgo nominal sensitivity curve which has been computed assuming a 10 W laser. 



the sensitivity at high frequencies. FIG. [T] displays and compares the progress achieved on the sensitivity curve since 
the beginning of the Virgo commissioning. At high frequencies (> 300 Hz), the sensitivity was limited by the shot 
noise and the laser frequency noise. At low frequencies (< 100 Hz) the longitudinal and angular controls of the mirrors 
are the main limitation; the electronic noise of the actuator and/or the sensor introduced by the control loops induce 
a displacement of the mirror which is well above the fundamental noise floor. 

2.2. C7 data set 

After the lock acquisition is completed, the interferometer state is set to science mode. The science mode duty cycle 
maintained by Virgo during C7 was 66%. The data was split into 53 data blocks, or segments, spanning over the 5 
days. A minimal duration of 10 minutes has been required for a segment. In order to define the start and end time 
of the data segments of good quality for the search for burst GW events we considered all the known instrumental 
effects that influence and spoil the output of the burst pipeline. For example, the last 10 seconds of each segment has 
been removed because many signals, including the GW strain channel, start oscillating and thereby cause the loss of 
lock. Another problem concerns the excitation of some thermal resonances of the last stage of the suspension (violin 
modes) when attempting to lock the interferometer. The decay time of these resonances has been determined, and 
for each segment we suppressed the beginning until the amplitude of the resonance having the longest decay time has 
been decreased by 90%. The reduction of the size of the data segments (hereafter called science segments) due to 
these data quality criteria amounts to 8%. In addition, some signals have been artificially injected into the GW strain 
channel during C7 by applying a force on the input mirror of one of the interferometer arm's Fabry-Perot cavity. 
These are the so-called hardware injections. The hardware injection signals are used in this analysis to test that the 
vetoes against external sources of noise that were implemented do not suppress any potentially real GW events (see 
Section r3.4.3p . The hardware injection periods were not considered for the GW burst search. Taking into account 
the application of the data quality criteria and the hardware injection periods, the duty cycle is finally reduced from 
66% to 55%. This corresponds to a total duration of 2.51 days. 
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3. DESCRIPTION OF THE BURST SEARCH PIPELINE 



3.1. Overview 



The GW burst pipeline that has been used in this search is composed of several parts that allow us to select and 
analyze all segments of Virgo data and eliminate, as much as possible, the artifacts in the detector output that could 
mimic a GW event. The core of the pipeline is the GW detection statistic, described in Section [3T3l that is applied on 
pre-selected segments of data of good quality. This results in lists of triggers. As we are dealing with non Gaussian 
and non stationary data, many of these triggers are due to instrumental effects. The next step consists of identifying 
all the sources of noise disturbances that pollute the detector output, and then define vetoes to be applied a posteriori 
on the trigger lists. This step in this pipeline is playing an important role as we cannot suppress these artifacts by 
requesting coincidence with another detector. We provide details on how the vetoes have been defined and which 
sources of noise they help in suppressing. 



3.2. Data selection 



Instrumental problems and environmental conditions have been identified that temporarily affect the detector 
sensitivity for a GW burst search. In a pre-selection step, we discard the periods where the detector operation is not 
optimal. To do so, we have defined a list of data quality fiags (DQ) that highlight the intervals of malfunctioning. 
Among the DQ fiags, some pertain to the saturation of the photodiodes and/or actuator signals that are used to 
keep the interferometer locked to its operating point. The analog electronics of the control loop that maintains the 
laser frequency noise below the requirements was sometimes also saturating. That then induces incorrect control 
signals and consequently a misbehaving GW strain amplitude. An excess of events was found when the GW strain 
h(t) reconstruction process was facing a problem, for instance when the calibration lines were buried in the noise. 
When this happens the h{t) reconstruction program sets a DQ flag stored in the data stream. Sometimes during C7 
a few frames or channels were missing due to some failure of the data acquisition system. That produces a hole in the 
continuous time series. The burst pipelines manage properly with the presence of holes in a segment, but to accurately 
compute the effective duty cycle we have defined a DQ flag pertaining to this issue. Finally, it has been found that 
when an aircraft is flying above the interferometer at a relatively low altitude acoustic and/or seismic noise couples 
into the dark fringe beam and induces a strong effect in h(t) that is detected by the burst pipelines. The band limited 
RMS of some acoustic probes (located in the end arm buildings) were monitored to detect the airplanes. All these 
DQ segments have been combined into a single list, taking into account possible overlaps. The dead time induced by 
the application of these DQ fiags on triggers list was only 0.8%, of the total duration which is acceptably low. 



3.3. Event trigger generation: Exponential Gaussian Correlator 



Different pipelines have been developed to search for GW bursts. Many of them [32, l33|, |3J, |35| rely on a common 
principle i.e., the detection of clusters of excess energy in a time-frequency representation of the data. The pipeline 
we use in this analysis follows the same idea; we refer to it as the Exponential Gaussian Correlator (EGG) as it been 
described in [36| . but here we give here a short summary. 

The time-frequency plane can be tiled using a lattice of sine Gaussian waveforms i.e., 

$(t)=exp|^--(^-j je^"/»*, (1) 

where the values of the central frequency /o and typical duration tq can be chosen in such way as to optimize the 
coverage of the plane. The idea is to consider the above waveforms as typical burst transients and to search for them 
using a matched filtering technique. For this reason, we designate ^(t) as a template waveform. The EGG computes 
the cross-correlation of the data with the templates, namely 

1 /•+"*(/)*•(/) .2../., 




where x{f) and $(/) are the Fourier transforms of the data and template, and >5'(/) is the two-sided power spectral 
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density of the noise. N is the normaUzation factor of the templates defined as 

For the present analysis, the power spectral density has been estimated over data segments of 600 second duration. 
To define the template lattice, the parameterization using the quality factor Qq = 27rTo/o instead of the duration tq 
is preferred. Using the algorithm of [s^l, we tile the parameter space (/o,Qo) with a minimal match of 99%. The 
parameter ranges 150 Hz < /q < 2 kHz and 2 < Qq ^ 16 correspond to the frequency band of best sensitivity for Virgo 
during the C7 run. T his genera tes 420 templates; the shortest and longest ones have 0.31 ms and 34 ms duration. 

The quantity p = yjl x |Cp is the signal-to-noise ratio (SNR), which we use as a detection statistic. It depends on 
the analysis time t^the template frequency /o, and quality factor Qo- It thus defines a three-dimensional representation 
map of the data [36] in which we search for a local excess as compared to typical noise fluctuations. 

We first apply a low threshold (SNR of 3.3) to the map. We form clusters with the surviving pixels of e nerg y |Cip 
with a two pass procedure. The pixels are grouped into clusters using the Hoshen-Kopelman algorithm [3^|. Once 
the cluster is formed, all pixels that have a SNR lower than 5 are removed. Then, all clusters overlapping in time or 
separated by less than 50 ms are grouped together once again. We produce the final list of triggers by requesting that 
the cluster SNR (\/2 x ^ |C'ip) is larger than 11.3. This choice for the threshold is a good compromise between not 
having too many triggers and not losing too much detection efficiency. 

The triggers' information is extracted from the clusters. The peak time, peak frequency and the trigger SNR are 
defined by the pixel that has the highest SNR. The duration and the frequency bandwidth of the trigger take into 
account all the pixels in the cluster. The timing resolution of EGG has been estimated for different kinds of waveforms 
[1^, |3^. For most signals used in the benchmark test, the timing resolution of EGG is smaller than 1 ms for SNR of 
10, but it can be as large as a few ms. 

3.4. Vetoes 

The search for GW burst events presented here was carried out on real data that are non-Gaussian and non- 
stationary. This prevents the use of theoretical estimates of the expected false alarm rate. On the contrary, the 
background rate must be extracted from the data itself. In the multi-detector context, the background can be 
estimated by c omp uting the rate of (purely accidental) coincidences between trigger lists that were time shifted with 
random delays [IJIl^l. We cannot resort to such procedure in the case of this study since we only have the Virgo G7 
data set at our disposal. There is a need to initially identify, in an thorough way, all the sources of noise that generate 
large SNR triggers that can mimic GW burst events; these types of noise events have to be suppressed before we try 
to estimate the background event number. 

We define a methodology that guaranties that we do not discard any true GW event as long as they do not coincide 
with a bad quality data period. Initially, there was the study of all the sources of noise that generate large SNR 
triggers, or more generally an excess of triggers, and the construction of vetoes to suppress these triggers. The study 
was performed on a subset of the G7 dataset (the playground segment) chosen arbitrarily before starting the GW 
burst search, in order to avoid bias. Our main concern in the veto development process was to minimize the dead 
time introduced by the data quality selection, while assuring a high level of noise trigger rejection. To perform the 
background studies and develop the veto strategy we used the longest G7 science segment; this segment was 14 hours 
in duration (22% of the whole dataset used for the GW burst search). 

Several transient detection algorithms have been applied to the playground dataset to identify and characterize the 
sources of noise that generate the high SNR triggers that mimic a GW burst event. Mean Filter searches for an excess 
in a moving average computed on whitened data [40]. Peak Gorrelator is a matched filter using Gaussian waveform 
templates [40]. The Wavelet Detection Filter is based on a multi- resolution time- frequency transform applied on 
whitened data [H]. Results of these pipelines have been compared to ensure that all sources of transient events have 
been discovered. These detection algorithms search for a short duration excess of energy, but using different methods. 
All of the trigger lists show a rather large excess of events compared to expectation from Gaussian noise; there was a 
definite excess of large SNR triggers. 

Two categories of large triggers have been identified; the first one corresponds to short glitches correlated to a glitch 
also present in an auxiliary channel. These auxiliary channels includes the available environmental monitoring and 
control loop signals (see Section [3.4. ip . The second category of glitches have been found to be due to the so-called 
"bursts of bursts" that we describe in Section [3.4.21 Specific vetoes against these two sorts of glitches have been 
developed as explained in Section r3. 4. 31 
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3.4-1- Origin of the short duration glitches 

To veto all periods of data during which instrumental or environmental problems occurred and generated a glitch in 
the GW channel data, one first needs to figure out which auxiliary channels are the most sensitive to these particular 
sources of noise. The by-eye scanning of the loudest events observed in the playground dataset already gave hints to 
potentially useful auxiliary channels. However we subsequently performed a systematic analysis using many channels 
recorded by Virgo and using the Mean Filter pipeline. 

We have found that three channels were particularly interesting for vetoing glitches in the GW strain amplitude 
channel h{t). It has been seen that many high SNR burst events are coincident in time with the realignment of the 
quadrant photodiodes by some stepping motors. These quadrant photodiodes, used for the interferometer automatic 
alignment system, are located at the output of the interferometer, on the same optical bench hosting the photodiodes 
used to detect the interferometer dark port signal. The stepping motors were generating very loud acoustic and 
mechanical noise that coupled into the dark fringe beam nearby. Each time a quadrant photodiode is moved there is a 
large increase in the signal measured by an accelerometer probe located on the detection bench. When the quadrant 
photodiodes were not re-aligned (horizontal and vertical positions remain constant) the RMS of the accelerometer 
probe signal remained low. The rate of the quadrant photodiodes' realignment was not constant over the run; initially 
the alignment occurred every 3 seconds, but this was then reduced to every 300 seconds at the end the C7 run after 
the discovery of the problem (the parameters of the process responsible for the quadrant photodiodes realignment 
have subsequently been changed to diminish this problem). This acoustic/mechanical noise generated glitches in the 
GW strain amplitude channel whose frequency content was around 550 Hz. Due to the fact that the excess of acoustic 
noise could last up to 1 second, which is long for a glitch finding algorithm, we used the RMS computed over 100 ms of 
data from an accelerometer located on the detection bench instead of the Mean Filter triggers to create a veto (called 
'Seismic' in the writing below). This allowed a more precise definition of the GPS time of the maximum of the excess 
of noise. This source of glitches has been found to be dominant in the playground dataset, but other (and rarer) 
types of glitches have also been found in the C7 data. Indeed, events in the Second Stage Frequency Stabilization 
System correction signal were observed in coincidence with h(t) glitches. Some remaining loud glitches discovered 
in the GW strain amplitude channel have been found to be coincident in time with a dip in the power of the beam 
reflected by the interferometer toward the laser as recorded by the so-called B2 photodiode. The light impinging on 
the B2 photodiode is sensitive to the power recycling cavity length change. The origin of the power dips have not 
been understood, but the effect on the GW strain amplitude was demonstrated and the safety of this veto using an 
optical auxiliary channel has been carefully studied, and is discussed below. To define an event by event veto for the 
two latest categories of glitches we used the Mean Filter triggers. These vetoes are hereafter referred to as 'SSFS' and 
'B2'. 



3-4-2- Bursts of bursts 

We noticed that the burst triggers were not uniformly distributed over time but seemed to show up in bunches 
lasting up to a few seconds. These events have been called Bursts of Bursts (BoB). FIG. [2] shows the whitened GW 
strain channel's amplitude and its spectrogram during a data period that contains a burst of bursts. 

The spectrogram shows that the event has a broadband frequency content contrary to most of the events caused 
by external disturbances. This broadband spectral signature provided clues that helped to identify the BoBs as local 
increases of the noise level due to an increase of the coupling factor between the frequency noise (which is only a factor 
2 lower than the shot noise at high frequency for C7) and the interferometer's dark port channel. The presence of a 
peak around T ~ 27 s in the correlogram of the burst events (see FIG. [3]) provided evidence that the residual angular 
motion of the mirrors could be playing a role. Indeed, 27 s corresponds to the period of a mechanical resonance 
between the two last stages of the Virgo suspension that, if excited, may induce mirror angular degrees of freedom 
motion (ilj . 

It turns out that the mirrors' alignment control system was not optimally working during C7 and the angular tilts 
of the mirrors induced an increase of the coupling of the laser frequency noise with the dark port amplitude. The 
BoB events represent a rather large fraction of the burst triggers in C7, as can be seen in FIG. [31 where the number of 
triggers in the correlogram fluctuates by up to ~ 50%. This is also visible by the height and width of the 37 mHz peak 
in the Fourier transform of the correlogram. The BoBs contribute partly to the highest SNR triggers, but mainly to 
the moderate SNR values. 

These non-stationary noise periods must be eliminated from the analysis since they generate an excess of noise 
triggers and they correspond to periods where the Virgo sensitivity is degraded compared to the norm. To identify 
those periods we used the band limited RMS of the dark port channel amplitude around llllHz, which is the 
frequency of the line injected to measure the laser frequency noise component in the dark port channel. This line 
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FIG. 2: An example of a Burst of Bursts (BoB) around t ~ 10 s. The upper plot shows the whitened interferometer dark port 
channel amplitude as a function of time, while the bottom plot is the corresponding spectrogram. 



„ 0.25 



T3 

c 
o 
o 
a> 

c 



> 




0.15 



0.05 



20 30 40 

Correlogram [s] 



60 




10" 

frequency [Hz] 



FIG. 3: Left: the correlogram of the triggers found by the EGG pipeline in the G7 dataset. A peak around T ~ 27s is present, 
with a replica at 54 s. A 0.6 Hz modulation is also visible, and this frequency corresponds to the pendulum mode of the last 
stage of the Virgo suspension. Right: the Fourier transform of the correlogram. The 37 mHz mechanical resonance between 
the two last stages of the Virgo suspension is clearly visible. 



allows one to measure the common mode noise level in the GW strain amplitude, and especially the laser frequency 
noise. Actually, measuring the height of the line allows one to monitor the variation of the common mode coupling 
factor [i^]. We computed the band limited RMS of the dark port channel amplitude around llllHz over 10 Hz 
(1106-1116 Hz), at a sampling rate of 20Hz. The high values of the band limited RMS correspond to periods where 
the non-stationary excess of noise is large. That allowed us to define the 'llllHz' veto. 
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3.4-3. Veto parameters tuning 

Once the channel and filter have been identified, one needs to tune the parameters of the veto, namely the threshold 
and the duration of the vetoed window around the time of the glitch. The parameters are adjusted to maximize the 
veto significance, keeping the dead time (i.e., percentage of vetoed science time) below a reasonable value. The 
significance is defined as the ratio of the number of triggers (presumed to be instrumental glitches) that the procedure 
vetoes by the square root of the number of triggers it would accidentally veto if there is no physical link between the 
glitches of the auxiliary channel and h{t) ^H]. The significance measures the excess of coincident events in standard 
deviation units. We used the GW strain amplitude triggers of the EGG pipeline to tunc the vetoes as we wanted an 
optimal background rejection for this analysis. The rate of accidental coincidence was measured by artificially time 
shifting the list of auxiliary glitches while the number of vetoed triggers is obtained for the zero-lag. We vary the 
threshold on the auxiliary channel output (SNR or RMS) and the duration of the veto window. For each value we 
compute the dead time, the efficiency (fraction of EGG triggers that are vetoed) and the use percentage of the veto 
(fraction of auxiliary channel triggers that veto an EGG trigger). A good veto must have a high use percentage [isj . 
The size of the veto window must be larger than the peak time difference in order to conservatively veto the region 
around the glitch, and can be as large as the total duration of the auxiliary channel trigger. Values from 100 ms up 
to 600 ms have been tested. 

FIG. [4] shows the number of coincident triggers between the GW strain amplitude and the auxiliary channel that 
has been used to define the Seismic veto, as a function of the time shift introduced in one of the lists. In this case, 
the auxiliary channel triggers are the RMS of the channel computed at 10 Hz. Different values of the threshold have 
been tested (from 150 up to 650 by steps of 50). The maximal value of the significance is obtained for a threshold of 
300. The threshold is confirmed by the fact that below this value the veto efficiency reaches a plateau. The choice of 
the duration of the veto window is based on a compromise between a high use percentage and a low dead time for the 
value of the threshold chosen. We checked that the optimal value, 500 ms, matches the width of the time difference 
between the auxiliary channel RMS and the EGG triggers. 
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FIG. 4: Number of coincident triggers between the GW strain amplitude and tlie auxiliary channel tliat lias been used to define 
the Seismic veto, as a function of the time shift introduced in the Seismic veto list. The different curves have been obtained for 
different threshold values (T) . S is the significance. The peaks correspond to the stepping motor that have been activated each 
6 seconds during the playground dataset. The different curves have been obtained for various values of the threshold applied 
on the auxiliary channel RMS. The maximal value of the significance is obtained for a threshold of T=300. 



The same procedure has been applied to tune the parameters of the 'SSFS' and 'B2' vetoes using the playground 
dataset. Some differences, compared to the Seismic case should be noted. For these two vetoes, built on the triggers 
generated by the Mean Filter glitch finding algorithm, wc used the information of the trigger duration to define 
the veto duration, in addition to a fixed minimal window size. For the llllHz veto, the duration of the window 
was defined by the time when the band RMS was above the threshold. The veto duration was then symmetrically 
enlarged. Given the nature of the BoB events (high EGG trigger multiplicity and rather noisy periods) we could not 
use the significance to optimize the threshold as it leads to setting the threshold at a very low value corresponding 
to an unacceptable large dead time. Instead, we set the threshold such that the dead time remained below 20%. 
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We checked that this choice assures that the veto was able to suppress a large fraction of the BoB events. The veto 
definitions and the chosen parameters are reported in TABLE |T1 



Veto 


Auxiliary channel quantity 


threshold tuning method 


threshold 


time window 


lUlHz 


Band RMS (computed at 20Hz) 


dead time 


12 


period above the threshold +/-300 ms 


Seismic 


RMS (computed at lOHz) 


significance 


300 


500 ms (fixed) 


SSFS 


SNR of the Mean Filter triggers 


significance 


6 


Mean Filter trigger duration - minimal value: 300ms 


B2 


SNR of the Mean Filter triggers 


significance 


10 


Mean Filter trigger duration - minimal value: 300ms 



TABLE I; Definition and parameters of the four event by event vetoes that have been developed and tuned for this analysis. 



3.4-4- Veto safety 

The safety of a veto is of fundamental importance as we do not want to inadvertently eliminate any real GW event. 
Vetoes based and dedicated to suppress coincident glitches in the GW strain amplitude channel must be safe with 
respect to real GW events. The DQ flags and vetoes which suppress deleterious periods of data can be potentially 
unsafe since the source of noise is independent from the effect of a real GW impinging upon the detector. On the 
contrary, we must be sure that a veto based on an auxiliary channel that has some glitches remains silent when a 
real GW event is passing through the detector. More precisely, environmental channels such as acoustic or seismic 
probes are expected to be safe, whereas vetoes constructed on optical or interferometer control channels may be 
unsafe, since a GW event will generate a change in photodiode signals that can be used within the feedback loop 
of some control systems. To test the safety of an auxiliary channel, one can examine the periods where hardware 
injections of fictitious GW signals are inserted into the interferometer. A deterministic force applied on the North 
Input (NI) mirror, induces a variation of the length of the North arm Fabry-Perot cavity which then mimics the effect 
of a real GW on the interferometer. During C7 there were two periods of hardware injections, and these included 
two different type of burst waveforms: 60 Sine Gaussian {f—920 Hz, Q=15, and /=460 Hz, (5=15), and 33 Gaussian 
(cr=lms). Both of these signals were injected with a SNR of 15 using a sensitivity curve taken just before the run 
started (the real SNR of the hardware injections is somehow different due to the sensitivity variation during the run). 
The critical point concerns the safety of the SSFS and B2 vetoes, which could in principle be unsafe. None of the 
hardware injection signals have been vetoed by the Seismic, SSFS and B2 vetoes. This establishes that our vetoes were 
safe. The anti-BoB veto llllHz which was built using directly the gravitational wave channel was, by construction, 
unsafe for GW burst event whose frequency content is around 1111 ± 5 Hz. We could not check that the veto was 
unsafe for these signals using the hardware injections as no signal with enough amplitude in this frequency range has 
been injected, but we should consider that the GW burst analysis presented here is insensitive for GW signals in the 
1106-1116 Hz frequency band. Besides, the safety of a veto assumes that its dead time remains small. This is the 
case for all the vetoes used in this analysis except the llllHz veto whose large dead time (16.1%) might suppress a 
real GW burst event. The anti-BoB veto llllHz suppressed 2 hardware injections, which is less than what we could 
foresee given the large dead time of this veto. 

4. OUTPUT OF THE GW BURST SEARCH PIPELINE 
4.1. First results 

We have generated the veto lists for the full C7 data set applying the tunings that have been defined using the 
playground segment as explained in Section [3.4.31 After applying the vetoes on the whole dataset, the remaining 
largest SNR triggers were studied in order to assess their compatibility with a GW burst signal. At this stage, we 
ended up with a single rather large SNR event that has been subsequently thoroughly studied. This very peculiar 
event has been observed outside the playground segment. Only one event of this type has been found in the full run, at 
least with such a SNR. It looks like a Sine Gaussian signal, as shown in FIG. [51 with a SNR of 69.7, for one of the Sine 
Gaussian templates (f=1466Hz, Q=15.5, a — (3/H/=3.4 ms); All environmental channels have been checked around 
this period using different glitch finding algorithms. Nothing suspect was found in any of the channels. However, the 
event was visible in the other phase of the demodulated dark port signal. We concluded that this event could not 
have been generated by a GW burst, but more likely could be due to an experimental problem, not occurring in the 
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chosen playground segment. Indeed, the modulation phase angle was tuned such that the effect of a GW crossing the 
interferometer was contained in one demodulated phase (ACp) of the dark port signal, while the other (ACq) should 
not be perturbed. This observation lead us to develop an event by event veto (called 'PQ') based on the ratio of 
the SNR of time coincident triggers in the two demodulated phase signals as proposed in the literature [H, |4^, . 
Despite the fact that such kind of events have not been found in the playground segment as a significant source of 
loud noise events, we decided to develop this veto and to apply it on the full data set. The other loudest high SNR 
events have been found compatible with BoB-like events. 
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FIG. 5: Time series in the dark port channel (left) and the whitened dark port channel (right) of the loudest event detected by 
the EGG pipeline in the search for GW bursts in the G7 data. This event was seen with a SNR of 69.7 at a frequency about 
1466 Hz. 



4.2. Necessity of an a posteriori veto 



The in-phase signal contains a priori the GW strain amplitude, provided that the demodulation phase is well tuned. 
An error in the demodulation phase induces a small coupling of the GW signal with the quadrature channel. However, 
the ratio of the GW energy seen in the two phases is expected to remain high. The ratio is expected to be proportional 
to ; where 54> is the error on the demodulation phase (45| . Unfortunately we do not know precisely the value 

of 54> during C7. A real GW event will be seen with a SNR in ACp much higher than in ACq. On the contrary, ACq 
will be sensitive to glitches in signals which are related to common mode noise. It may happen that some source of 
noise can affect both quadrature signals with similar strength. This is especially the case for a dust particle crossing 
the laser beam before the output dark port photodiode where the beam is especially small in width. A key point 
of such a veto is to assure that no real GW events would be suppressed, and therefore one would need to develop a 
veto with a rather good security factor. To do so, we used the hardware injections to verify safety and develop the 
characteristics of this veto. 

EGG triggers have been generated for the dark port ACp and ACq channels and all events with SNR above 5 
were considered. We then required time coincidence between the two sets of triggers using a window of ±10 ms. No 
coincidence in frequency has been required. The hardware injection signals were detected in the ACp triggers with an 
efficiency of 98%. Unfortunately during C7, the injected SNR was not, on average, very large and consequently the 
signal energy in the ACq channel could not be expected to be very high (below the pipeline threshold at SNR=5). 
Among the 91 detected hardware injections only 5 had a coincident trigger in the ACq phase channel within 10 ms; 
this tends to indicate that a large fraction of the GW energy would be contained in the ACp channel as expected; this 
is displayed in FIG. El The ratio between the SNR detected in the two demodulation phase channels in coincidence is 

^ SNRacp 
" SNRAc^i 
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ACp SNR 

FIG. 6: SNR of the coincident EGG triggers in the two quadrature demodulated channels of the output dark port signal: ACp 
(in phase) and AGq (in quadrature). The time coincidence window of the triggers seen in the two channels is 10 ms. The 
Seismic, SSFS, B2 and llllHz vetoes have already been applied. The peculiar event appears very isolated with respect to the 
other triggers with its ACq SNR of 141. The hardware injection signals, which were seen in both of the two demodulation 
phase channels, are indicated by the circles for burst-like signals and by the squares for the compact binary coalescence (CBG) 
signals. 



The very low event statistics (5 events detected simultaneously in ACp and ACq) encourages us to stay conservative 
when defining the PQ veto parameters. We checked that the use of a ±10 ms window limits the accidental coincidence 
event rate to 0.2 over the period of the hardware injections (assuming a Poisson trigger rate). This excludes the fact 
that we have more than one accidental association with an ACq trigger for the fix burst hardware injections as shown 
in FIG. El In this figure, the strange and unique high SNR event is totally isolated from the rest of the events. This 
event has k ~ 0.42 and a SNR in the ACq channel of 141; these facts totally exclude the possibility that this event 
has been generated by a real GW. 

The goal of the PQ veto is to suppress high SNR events in the ACq channel that induce a transient in the dark 
port channel amplitude. That is why it is reasonable to consider only large SNR ACq triggers. A threshold is applied 
on the SNR of the ACq triggers {SNRacq)- The other veto parameter threshold pertains to the ratio k. We did 
not try to optimize these two parameters, but rather we keep a conservative attitude with respect to the hardware 
injection signals' position in the two-dimensional plane that shows the SNR of both quadratures for the coincident 
triggers (see FIG-H]). The veto list segments have been defined as the periods during which SNRacq > 8 and k < 1. 
All burst hardware injection signals pass these conditions with a safety factor since they are all below the k — 2 line. 
The starting and ending times of the veto segments are given by the ACq trigger times. The dead time of the PQ 
veto is 0.07%, which is very acceptable. 

4.3. Final results 

TABLE im gives the dead time of the five event- by-event vetoes which have been used in this analysis. The total 
dead time amounts to 21.2% taking into account overlaping veto segments. This is a rather high value, mainly because 
of the presence of the non-stationary excess of noise due the looseness of the mirror's angular degree of freedom control 
during C7. 

FIG. [7| shows the triggers' SNR distribution obtained from the full C7 data set, after having applied consecutively 
the five vetoes (note that the highest SNR event is off of the scale of the plot). The four anti-glitch vetoes (Seismic, 
SSFS, B2 and PQ) suppress 10.6% of the total number of triggers and play an important role for the large SNR 
(>10) events since 51.3% of them are eliminated. The anti-BoB veto is critically important, and its efficiency on EGG 
triggers is impressive; 45.6% of the remaining triggers are eliminated, and 97.6% of the remaining loudest ones (SNR 
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Veto name 


Seismic 


SSFS 


B2 


PQ 


llllHz 


all vetoes 


dead time 


3.75% 


2.6% 


0.03% 


0.07% 


16.1% 


21.2% 



TABLE II; Dead time of the diflerent vetoes that have been found useful. The values are computed for the full C7 dataset. 
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FIG. 7: SNR distribution of the triggers obtained on the GW strain amplitude channel by EGG after applying all vetoes on 
the full C7 data set. 



> 10). FIG. [H] shows the SNR of the remaining triggers as a function of their frequency. The application of all these 
vetoes allows us to reduce, by a factor 2, the number of events for all EGG triggers above SNR=5. 98.8% of the EGG 
events above SNR>10 were eliminated. 21 triggers above SNR=10 remain at the end. The loudest has a SNR of 
11.26, and although it was not vetoed, it seems to be due to a short period non stationary Virgo noise increase. 

5. GW BURST SEARCH ANALYSIS SENSITIVITY 
5.1. Generation of the injected waveforms 

The GW burst search presented in this paper is intended to be efficient for all waveforms as long as they have energy 
concentrated in time and frequency in the range 150-2000 Hz. However, since the Virgo sensitivity is not constant over 
this frequency band, the detection efficiency for a given signal strength will depend on the frequency of the signal. In 
order to assess the sensitivity of the detector to dissimilar kinds of signals at different frequencies, we used simulated 
waveforms of diverse types. These were added to the Virgo GW strain amplitude (software injection) and analyzed 
by the EGG pipeline, similar to how the G7 triggers were produced. These injected signals were randomly distributed 
in time during the G7 run, assuring a minimal separation of 60 s [13] • 

The waveforms investigated in this analysis were of two types. One group consisted of ad-hoc waveforms, such 
as Sine Gaussians and Gaussians, whose detection efficiency is expected to depend mainly on the central frequency, 
bandwidth and duration of the waveform. On the other hand, to test the efficiency of this analysis to detect BBH 
merger and ring-down GW signals we used waveforms provided by numerical relativity simulations (48j . We also used 
a signal pertaining to a star core collapse, referred to as A1B2G1, from the supernova simulation catalog available in 

m- 

The amplitude of the GW associated with an astrophysical source depends on the distance of the source from the 
Earth. To quantify the strength of a GW burst, we used the root sum squared of the strain amplitude at the Earth 
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FIG. 8: SNR as a function of the frequency of the triggers obtained on the GW strain amplitude channel by EGG after applying 
all vetoes on the full C7 data set. The distribution around 2 kHz is due to the templates' placement; some Sine Gaussian 
templates have a central frequency which can slightly exceed 2kHz. 



without folding in the detector antenna pattern 



hrss = ^ J {\h^it)\^ + \h+{t)\^)dt, (4) 

where and hx are the two GW polarizations that we express below for the different cases relevant to our analysis. 
We assume in Eq. ^ that the source is optimally oriented with respect to the detector. 



5.1.1. Sine Gaussian and Gaussian waveforms 



For the Sine Gaussian and Gaussian signals, the detection efficiency is expected to depend mainly on the central 
frequency, bandwidth and duration of the waveform. For these signals the strain sensitivity h{t) of the detector can 
be written as 

/i(t) = + (5) 

where and are the antenna response functions characterized by the source position (0, 0) and polarization 
angle \1/ relative to the detector Specifically, for Sine Gaussian injections we used circular polarized waveforms 
described by 

hAt)^h,, e-2(-/ot/Qo)= sin(2^/ot), 

h+{t) = /lo e-2(-/o*/Qo)^ cos(2^/ot). (6) 

The typical duration of a Sine Gaussian signal depends on the parameter Qo and the central frequency /o as ^^^^ . 
Several values of /o ranging from 100 Hz up to 1797 Hz and (3o=3, 8.9 and 9 were used. Their narrow-band feature 
allows the testing of the pipeline performance in a given frequency region. 
For Gaussian injections, we used linearly polarized waveform given by 



h+{t) = hoe-'"/^\ 
hx = 0, 



(7) 
(8) 
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where r is the width of the signal. Gaussian waveforms are of interest in the GW burst search because many predicted 
signals associated with core collapse GW emission just after the bounce have large peak structures. We chose the 
values of r to lie between 0.1 and 6 ms. In the case of linearly polarized Gaussians, the angles 0, 4> have been randomly 
chosen. For the circularly polarized Sine Gaussians the randomization has been done over the angles 9, cf) and ^E". 



5.1.2. Astrophysical waveforms 

For the core collapse signals, we have chosen only one waveform produced by a core collapse 3D simulation 
conducted with numerical general relativistic techniques in order to estimate, with a physical model, the detection 
efficiency for a realistic supernova GW signal. This signal corresponds to the case of a stiff equation of state, a 
small initial differential rotation, and a moderate rotational kinetic energy [49|; this leads to a waveform with a 
large negative peak followed by a ring-down phase (regular collapse). More up-to-date simulations [5l|, [sll using 
hydrodynamical models with realistic nuclear equations of state tend to confirm the general features of the regular 
core collapse. In this present study we are not interested in the details of the waveform. The waveform is linearly 
polarized and the simulations have been done similarly to the Gaussian signals. 

With the recent breakthroughs in the field of numerical relativity (NR) , many groups have been able to simulate 
the evolution of binary BHs through three stages: inspiral, merger and ring-down ^IS, i4S. i53j. These waveforms have 
been shown to be consistent with the already existing analytical results using post-Newtonian and BH perturbation 
theories f &l . [ssl . [sgI . Ist} for the inspiral and ring-down phases respectively. Furthermore, despite the differences in the 
numerical methods employed, gauges chosen, and methods adopted to evolve the systems, qualitatively the results 
from all groups show good agreement All these provide a motivation to use these results from a data analysis 
perspective. Many such attempts have already been made to search the three phases of the binary black hole evolution 
using NR outcomes [13, Ho] ■ We apply these results from the point of view of a GW burst data analysis and use 
the NR waveforms to assess the efficiency of the EGG burst pipeline to detect the BBH merger. For this application 
we used the numerical relativity (NR) simulations of the non-spinning, equal mass BBHs provided by the Goddard 
Space Flight Center group [i^. The injected waveforms were GW strains from the leading Z = 2, m = 2 spin- weighted 
spherical harmonics from the simulations. These waveforms had (approximately) 3 cycles of inspiral, followed by 
merger, and then a ring-down. 

Unlike the ad- hoc waveforms or the linearly polarized case of SN waveforms, one has to include the effect of the 
orbital inclination angle l for the BBH merger waveforms. Eq. ^ should be rewritten in order to include the effect 
of orbital inclination as 

h{t) = F+i0, 4>, ^) A+{l) h+{t- A - 0) + Fx (0, 0, ^) Ay, {i) hy {t- i = 0), (9) 

where and Fx are the usual antenna pattern functions, A+ = —^(1 + cos^ l), A^ — — cost and /i+(i; t = 0) and 
/ix(i; 1 = 0) are the polarizations for the situation of zero inclination angle t. In this case, the simulated waveforms 
are generated by randomizing all four angles involved. The NR waveforms for the BBH mergers in the total mass 
range 5Mq-15OM0 were injected into the C7 data in order to calibrate the detection efficiency of the pipeline to these 
types of mergers. The range of masses chosen for the injection was decided in accordance with the sensitivity of the 
C7 run. 



5.2. Detection efficiency of tlie GW burst search with the software injections 

We applied the same data quality criteria (DQ and vetoes) on the trigger list containing the waveform injections. 
The detection efficiency is defined as the fraction of injected signals reconstructed with a SNR larger than 11.3 
corresponding to the loudest event found in the C7 data triggers. As the veto dead time is rather large, the vetoed 
periods have not been considered when computing the detection efficiency. FIGs. El [10] and [TT] show the detection 
efficiency as a function of hrss for each family of waveforms. The error on the efficiency data points takes into 
account only the statistical error. To estimate the values of hrss at efficiency 50% and 90% the four parameters of an 
asymmetric sigmoid function have been fitted to the data points. The efficiency e(/irss) is defined as follows, 

mid :„ I -umid 



where e^ax is the maximal efficiency obtained for strong signals (it should tend to unity), r™ is the ratio hrss/h 
where ft-™/ is the strain amplitude at half height i.e., we have e(ft™/) = g™°^/2. a and /? are respectively the slope 
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FIG. 9: Detection efficiency of the search analysis as a function of the signal strength hrss at Earth for the Gaussian waveforms 
(with r specified in the legend) and one core collapse simulated waveform (A1B2G1). The efficiencies were computed for sources 
at random sky locations. The error bars take into account only the statistical error on the detection efficiency. 



and the asymmetry of the sigmoid function. The parameters were estimated by minimizing a function, e™'*^ is 
constrained to remain smaller than unity. For signals for which e™'*^ does not reach unity, we checked that this is due 
to a fraction of signals with an injected SNR smaller than the threshold of 11.3. This is especially the case for the 
Gaussian and the A1B2G1 core collapse simulated waveforms, as expected. Indeed, as the frequency content of the 
Gaussian signals is maximal at low frequency, for a fixed hrss and a random sky position, we expect a non negligible 
fraction of low SNR Gaussian signals. 

As expected for the Gaussian signals, the performance decreases as the width r of the Gaussian increases. For r 
greater than 4 ms, the detection efficiency is already quite poor, as shown in FIG. [HI Note that the results obtained for 
the supernova signal A1B2G1 are very similar to the Gaussian with t=0.25 ms (see TABLE irv|) : this is not surprising 
as the A1B2G1 signal has a large amplitude peak whose width is approximately 0.35 ms. FIG. [TOl shows the detection 
efficiency of the analysis obtained with the Sine Gaussian waveforms; the performance follows the trend of the C7 
Virgo noise floor shown in FIG. [TJ which is expected for a signal well localized in frequency, such as a Sine Gaussian. 

The performance obtained for the BBH merger waveforms for a total mass M from 5 Mq up to 150 M© (displayed 
in FIG. [TT|) were not as good as for the Sine Gaussian signals. The best h^^'^" was achieved for M=20 Mq with 
h^^f' ~ 2.2 X 10~^°/\/ Hz. That was at least two times worse than the best Sine Gaussian result. The merger 
frequency can be approximated by the empirical formula [T^jHIIj deduced from FIG. 8 in [6l|, 

_ 15kHz 

/merger — , (11) 

where M is the mass (in units of solar mass) of the flnal BH formed after merger. M is smaller than the sum of the 
initial BH masses. For M=20 AIq this corresponds to 750 Hz, which is well inside the frequency range of the search. 
On the contrary, for masses above 100 Mq, /merger is smaller than 150 Hz, which is the lower frequency boundary of 
this search. The EGG pipeline for this kind of waveform is still able to catch the ring-down part of the signal, which 
has an appreciable contribution to the SNR, while the performance for a Sine Gaussian at the /merger frequency would 
hardly be detectable. For the low mass signal the drop of performance is explained by the fact that the frequency at 
the merger exceeds the upper frequency boundary of this search, and moreover, the equivalent Q is expected to be 
larger than 16 (the upper value chosen in this search). 

The value of hrss corresponding to 50% and 90% of efficiency (/i^sf" and /i^s?) have been determined using the 
fitted sigmoid functions and are given in TABLES lllli IIVI and IVl for Si ne G aussians, Gaussians, Supernova, and BBH 
merger signals respectively. The best value of /i^s?! ^-'-'^ ^ 10~^°/V Hz, was obtained with the Sine Gaussian of 
frequency /o=361 Hz and Qq = 9; this signal is 15 times higher than the noise floor at this frequency. Note that bad 
performance obtained for the 100 Hz Sine Gaussian is expected since the lower frequency range of the search was 150 
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Sine Gaussians Q=3 




Sine Gaussians Q ~ 9 




FIG. 10: Detection efficiency of tlie search analysis as a function of the signal strength hrss at Earth for Sine Gaussian waveforms 
with different central frequencies /o and for two different values of the Qo factor: left Qo = 3, right Qo ~ 9- The efficiencies 
were computed for sources at random sky locations. The error bars take into account only the statistical error on the detection 
efficiency. 



Hz. Moreover the C7 sensitivity is quite poor below 150 Hz. 

The quality of the fitting procedure has been controlled by checking that the of the fitting function remains close 
to unity within at maximum a factor 3. Increasing the error on the efficiency by up to a factor 3 does not change the 
result on h^ss to within 1%. The statistical error on h^^^° and /i^^f" have been computed using the covariance matrix 
of the efficiency curves; for the Sine Gaussian and BBH waveforms the error on /i^^f" was below 1.3%. Note that for 
the errors were higher, between 6% and 13%. For the Gaussian and core collapse waveforms the uncertainty 
is on average larger, indicating that the quality of the fitting was worse. The errors are between 1% and 4% for all 
the Gaussians, except for the Gaussian with r=4 ms for which the error on reaches 11%. However, the biggest 
uncertainty on the hj^^f^'^" values comes from the error on the GW strain ampliude calibration, which is of 40%. 
Including the systematic error due to the calibration uncertainty leads to an error on h^^f^'^" between 40% and 43%. 



6. SEARCH RESULTS 



6.1. Event rate and GW strength upper limits 

From Section 0] we can assert that we have not observed any GW event with a SNR above 11.3, corresponding to 
the loudest event SNR. Knowing the detection efficiency as a function of the signal strength for different kinds of 
signals, we can now derive a 90% confidence level (CL) upper limit on the event rate that depends on the strength of 
the signal. Using the loudest event statistic as described in [b^I, the 90% CL exclusion on the event rate of a signal 
detected with an efficiency of e{hrss) is given by, 

^,0%(,^^^)^ZM1_^, (12) 
Te[hrss) 

where T ~ 2.1 days is the effective observation time for the C7 run (veto times have been subtracted from the 
total analyzed time). The results are given in FIG. [12] for the Sine Gaussian waveforms (Q ~ 9) and the Gaussian 
waveforms. 

As expected for strong signals, the exclusion rate reaches the asymptotic value of ~ 1.1 events per day 

whatever the signal type. We can exclude at 90% confidence level a GW event rate TZgQ%, assuming 100% detection 
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FIG. 11: Detection efficiencies from tfie searcli anaiysis as a function of the signal strength hrss for the BBH merger and ring- 
down waveforms for different BH total masses (but where the two individual component masses are equal). The efficiencies were 
computed for random polarization, inclination angles, and sky positions. The error bars take into account only the statistical 
error on the detection efficiency. 
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TABLE III: hrss corresponding to 50% and 90% efficiency in units of 10~^° strain /V Hz obtained for the ad-hoc Sine Gaussian 
waveforms. The efficiencies were averaged over the full C7 data set for random sky positions. The symbol - means that no 
simulated waveforms were available. The symbol (*) indicates that we used simulated waveform with 0=9 instead of 8.9. 



efliciency, of 1.1 events per day for T ~ 2.1 days of observation. On the other hand, depending on the characteristics 
of the signals the exclusion rate, when the signal becomes weak, is rather different. The vertical asymptotic value in 
FIG. [12] gives an indication of how the sensitivity of the search depends on a given frequency for the Sine Gaussian, 
or on the width of the Gaussian waveform. The exclusion results for Sine Gaussian waveforms with Q=3 (from 153 
Hz up to 1053 Hz only) are very similar to the Sine Gaussian results with Q=9. 
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r (ms) 50% efficiency 90% efficiency 
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12.4 



TABLE IV: hrsB corresponding to 50% and 90% efficiency in units of 10 strain /\/ Hz for the Gaussian waveforms of different 
widths and one core collapse simulated waveform (A1B2G1). The efficiencies were averaged over random sky position. 
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13.2 


32.4 
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24.5 


74.4 



TABLE V: hrss corresponding to 50% and 90% efficiency in units of 10 "^^ strain /V Hz for the BBH numerical relativity 
waveforms. The total mass of the system varies from 5 Mq up to 150 Mq. 

6.2. Astrophysical interpretation 

In the previous section we determined the sensitivity of the C7 data analysis to different potential GW burst sources 
in terms of hrss- This can be translated into an astrophysical estimate on the amount of GW energy emitted by a 
source or the distance reach, depending on the type of sources we consider. For simplicity we have assumed in our 
ad-hoc waveform simulations that the emission is isotropic, although this may not be the case for a real astrophysical 
signal. In the case of a Sine Gaussian waveform with a central frequency /o and quality factor Q S> 1, one obtains 
for the GW energy emitted by the source [gI, [6^ 

2 3 

i5GW=^^(27r/of /^L, (13) 

where r is the (non-cosmological) distance to the source. Similarly, for a Gaussian waveform, the emitted GW energy 
is 

t'^ c"^ 1 

^GW - /iL: (14) 

where r is the width of the Gaussian. 

The above expressions can be used to derive an astrophysical sensitivity with the present C7 data analysis in 
different ways. For sources at a fixed distance, an estimate on the emitted GW energy can be obtained by using 
the smallest value of /irsf"- Conversely, if we know the energy emitted in GWs, this information can be used to 
infer the maximum distance up to which such an event could have been observed by Virgo with an efficiency larger 
than 50%. In what follows we discuss these various possibilities, quantifying the detectability of supernova events 
and BBH merger events where the nature of waveforms and the energy of emission are known through the numerical 
simulations. 
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FIG. 12: Event rate as a function of the hrss exclusion plots at 90% confidence level obtained in the Virgo C7 data. The left 
plot corresponds to GW burst events that can be modeled by a Sine Gaussian with a Q of 8.9 or 9 for frequencies, /, covering 
the most promising part of the C7 frequency range. The right plot corresponds to GW burst events that can be modeled by a 
peak of different widths, r, from 0.1 ms up to 4 ms. 



6.2.1. Detectability of an arbitrary burst GW event close to the galactic center 

Considering the Sine Gaussian and Gaussian waveforms that give the best results in terms of h^^f' , we can compare 
the estimation of the maximal GW energy emitted by a hypothetical GW burst source in the vicinity of the galactic 
center, assuming a distance of 10 kpc. With the 361 Hz Sine Gaussian we obtained /i^o% of 1.05 x IQ-^o/ThJ, which 
translates into an energy of 3 x IO^^Mq (? . For a Gaussian waveform, the most favorable case corresponds to 
r = 0.25 ms and /i^s? — l-^S x 10^^"/\/Hz; the associated energy is ~ 2.8 x lO^'^M© c? . Hence the most optimistic 
limit on the energy differs almost by an order of magnitude between the Gaussian and Sine Gaussian waveforms, while 
the detection performance in terms of h^ss is only two times worse for a Gaussian with respect to a Sine Gaussian. 



6.2.2. Detectability of supernova events 



We now consider the predictions on the GW energy produced by the collapse of a stellar iron core to form a neutron 
star. It is assumed that the progenitor star rotates fast enough to generate sufficient deformation. We do not consider 
here the GWs generated after the core collapse by the development of an oscillation of the proto-neutron star lasting 
a few hundred of milliseconds [l0| since we did not expand the search for corresponding high Q oscillating signals. 
A few hundred of milliseconds at a frequency around 700 Hz corresponds to a Q of a few hundred. Using the core 
collapse waveform A1B2G1 injected into the data to determine the efficiency for that particular signal, we can derive 
the distance reach corresponding to 50% of efficiency as we do for the BBH waveforms; the distance is 0.15 kpc. 
However, this waveform might not be representative of the overall core collapse bounce waveforms. Indeed, depending 
on the progenitor model parameters (mass, rotation rate, equation of state) the GW energy prediction can vary by 
several orders of magnitude [52|: 3 x 10~^^ Mq < Eqw < 7 x 10~^Mq with peak frequency which can be as high 
as 800 Hz. However, these recent simulations (s^l seem to indicate that the waveforms are rather generic regardless of 
the dynamic of the collapse (pressure dominated or centrifugal bounce). Even in the case of centrifugal bounce, the 
waveforms exhibit a positive pre-bounce followed by a negative peak, but the typical frequency is lower than in the 
case of a uniform differential rotation for which the typical frequency is around 700 Hz. For instance, if one considers 
the model with a 20 Mq progenitor mass, a moderate differential rotation value (A=1000 km) and a moderate initial 
rotation rate (/3 = 0.5%), Eqw — 2.32 x 10^^ Mq and the characteristic frequency is 715 Hz using the equation 
of state provided by Shen [65;]. We should not specifically assume that supernova GW signals are monochromatic 
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waveforms like Sine Gaussians. They do not look like Gaussians either (they are more localized in frequency than 
Gaussian). Taking the h^°f' value for the r = 0.5 ms Gaussian (that roughly corresponds to the width of the peak 
of the considered waveform), we obtain a maximal distance range of 0.18 kpc. Using the Sine Gaussian expression, 
and using the /i^s? ^'^^ * central frequency 715 Hz we obtain a distance of 0.12 kpc. For the A1B2G1 waveform 
used in this paper, and considering a radiated energy of Eqw — 2.13 x 10~^Moc^ computed using Eq. (5) of 0], 
we obtain a distance of 0.12 kpc whatever the assumption made on the waveform (Gaussian of t = .35 ms or a Sine 
Gaussian of /o=650 Hz). These numbers indicate the order of magnitude of the maximal detection distance we have 
reached during Virgo's C7 run. This achieved distance is somewhat small, and we expect to gain at least one order 
of magnitude with Virgo data acquired during more recent observation runs. 

6.2.3. Distance reach for binary black hole mergers 

From the detection efficiency results obtained with the BBH waveforms in Section 15. 2i one can derive a distance 
at which EGG is able to detect a signal with 50% or 90% efficiency, respectively d^^'^" and These values have 

been extracted from the results of the fits of the detection efficiency as a function of the distances for which the BBH 
sources have been simulated. An asymmetric sigmoid function has been used and the errors on (f'^'^" and have 
been estimated as they have for /i^^f" and h^^^° . FIG. [T3l shows the and detection distances as a function 
of the total mass of the system. gives an indication of the average distance range up to which this search could 

have been able to detect a BBH source. With the present choice of threshold on the SNR, our search is sensitive to 
BBH mergers in the mass range 20 M0-15O Mq up to a distance of 1-3 Mpc (0.5-lMpc) for 50% (90%) detection 
efficiency. The distance reach is maximum for a BBH of total mass 80-90 M©. In the C7 sensitivity band, this mass 
range corresponds to signals where the merger part is dominant. Besides, for a mass higher than 100 Mq, the merger 
frequency (given by Eq. pT|) ) is below the frequency lower bound studied in the analysis (150 Hz). This fact adds to 
the degraded C7 run sensitivity below 150 Hz, and explains why the distance reach for the high mass waveforms is 
decreasing. 

Merger frequency [Hz] 

300 150 100 
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B 1.5 
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'0 50 100 150 

Black hole binary total mass [M J 

FIG. 13: Detection distance of the black hole binary source with efficiency of 90% and 50% as a function of the total mass of 
the system. The curves with circles and squares have been obtained by estimating the EGC pipeline detection efficiency to 
numerical relativity black hole binary waveforms. The curve with diamonds was derived from the results obtained with Sine 
Gaussian waveforms assuming that the waveform when the two black holes merge resembles a Sine Gaussian. The error bars 
take into account only the statistical error on the detection efficiency. 

This bound on the distance can, in principle, be translated into an upper limit on the event rate, assuming a 
Gaussian distribution of BBH masses around 80 Mq. This analysis would require combining a model for the galaxies, 
their blue light luminosity distribution as a function of distance and mass, along with the antenna pattern of the 
detector; this has been proposed in ^66j| for the LIGO detectors. But given the limited volume of the universe that 



-e-NRe = 50% 
-B-NRe = 90% 
■❖-SGe = 90% 
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the present search is capable of sampUng, we postpone this more complete analysis until we analyze the most recent 
Virgo data with its improved sensitivity. Moreover, the results have been obtained only for non-spinning equal mass 
binary systems. 

An alternative way of deriving the distance reach for BBH mergers is to assume that the merger and ring-down 
phases of the BBH waveforms resemble Sine Gaussians and use the results obtained from the Sine Gaussian waveforms 
to derive the detection distance. In order to use Sine Gaussians to model the BBH mergers, we approximate the NR 
waveforms to be Sine Gaussians with peak frequency /o equal to the frequency at the merger (/merger) of the NR 
waveforms using Eq. (jlip following [24| . It should be noted that this definition of /merger is only approximate. 
Besides, it is not obvious how to assign a precise value of Q to the Sine Gaussian that approximates a BBH waveform. 
However, taking into account the typical duration r of the "merger" part of the waveform and using the relation, 
exact for a Sine Gaussian, r = (3/(27r/o), one could deduce that the Q is between 10 and 20 for the high mass BBH 
and as large as 40 for the low masses. This is larger than the boundary of Sine Gaussian templates bank used in this 
analysis. Hence, we may suspect that for the low mass (< 2OM0) the EGG pipeline is not optimally reconstructing 
the SNR of the NR waveforms. However, it is straightforward to show that the quantity which is used in the 

derivation, does not depend on Q for Sine Gaussians as long as Q is sufficiently large fTABLE IIIII confirms this). For 
all these reasons, in the discussion that follows, we consider a Sine Gaussian with Q ~ 9. In order to determine the 
distance reach one should make an assumption on the total energy emitted by the source and on the source emission 
pattern. According to NR simulations f6l|, ~ 3.5% of the total rest mass energy would be radiated away as GWs. 
Following [24], we can use Eq. and the values of /i^ss in TABLE Hill to derive the distances up to which such BBH 
merger events could be observed with an efficiency of 50% or 90%. The results for the distance reach at an efficiency 
of 90% are presented in FIG. [T3] together with the distance reach obtained using NR waveforms at the same detection 
efficiency. At intermediate masses and frequencies the distance reach obtained when BBH mergers are modeled with 
SG waveforms is about 60% larger than the one obtained when more realistic NR waveforms are used. A detailed 
comparison of the two waveforms shows that this difference is partly due to the broader spectral distribution of the 
NR waveform for which part of the signal energy falls out of the detector and analysis bandwidth. The remaining 
difference is due to the EGG pipeline, which by construction matches more efficiently with the SG waveforms. 

7. CONCLUSION 

In this paper, we have reported on a search for short duration GW bursts in the Virgo G7 data from within 
the 150Hz - 2 kHz frequency band; this corresponds to the best sensitivity achieved by Virgo during the G7 run. 
This search, carried out using the output of only one GW detector, required first the understanding and then the 
possibility of vetoing all identified sources of noise. Special care has been taken to assure ourselves that the procedure 
does not suppress a GW event, and that the dead time remains small. Unfortunately, the non stationary data features 
prevent a search during vetoed time periods, which amounts to an accumulated dead time of 20.2 %. Though this 
rather large dead time limits the astrophysical interest of the search, it has been possible to calculate a sensitivity 
for various GW burst sources. The best sensitivity that has been reached in terms of the square root of the strain 
amplitude at 50% of efficiency, /i^sf", was 1.1 xlO^^" /\/Hz at 361 Hz, a frequency region where Virgo's sensitivity 
was optimal. If one considers the whole set of waveforms studied in this paper, the sensitivity lies between 10"^*^ and 
IO^^^/VHz. However, one should note that the error on the sensitivity can be as large as 40% due to uncertainties 
in the calibration of the strain amplitude. Those numbers can be compared to GW burst searches carried out with 
LIGO detectors during the S2 data taking in 2003 [21}. Indeed, the G7 Virgo sensitivity is comparable to that of 
4-kni LIGO detectors during their S2 run; Virgo during G7 was a little bit more sensitive above 500 Hz and a bit less 
below. Not surprisingly, the search sensitivity obtained with G7 Virgo data is similar to LIGO S2 results for the same 
kind of waveforms: having values between 10"^" and 10~^^ / VHz. 

Besides the search of un-modeled GW burst waveforms, we have shown using recently computed numerical relativity 
waveforms that burst pipelines can be used efficiently to search for BBH merger GW emission. In the case of non- 
spinning equal mass system, a distance reach of 2.9 Mpc for 80 Mq total mass has been inferred from the G7 data 
assuming an efficiency of 50%. This mass range corresponds to when the BH merger part of the signal dominates over 
the inspiral part. That also corresponds to a frequency region where ground base interferometers achieve their best 
sensitivity. It is interesting to note that the sensitivity in the merger dominated frequency range is up by roughly a 
factor of 10 in the most recent Virgo data (acquired in the later half of 2007). This would correspond to a distance 
reach up to 30 Mpc (12 Mpc) at 50% (90%) efficiency. This distance reach improvement is not enough to tackle 
regions where the number of galaxies becomes interesting given the expected BBH event rate. Further upgrades, 
like Advanced Virgo, may be needed to have a reasonable probability of detection. Despite that, and given the 
improvement of the Virgo sensitivity in the low frequency region, we will expand the EGG pipeline parameter space 
to explore more efficiently from this region. Though matched filtering with a bank of templates, including inspiral, 
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merger and ring-down phases, might have a better detection efficiency for BBH mergers, it is important to complement 
the search with other techniques not based on tailored matched filtering but using generic templates, such as Sine 
Gaussians. This approach could provide a method that would be more robust at finding increasingly complicated 
signals that depend on a larger number of astrophysical parameters from the sources, such as a black hole's spin in 
the case of a black hole - neutron star system. 

Finally, one should add that the recent Virgo data taking (VSRl) will be analyzed in coincidence with LIGO S5 
data. That is especially interesting in the high frequency region (> 600 Hz) where the sensitivity of the four detectors 
are all limited by shot noise. This network of four non-aligned interferometers with comparable sensitivity should be 
able to detect GW signal emitted by a supernova anywhere in the galaxy. 
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